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ABSTRACT 

Context. We present angle-dependent, broad-band intensity spectra from accretion disks around black holes of 10 M Q . In our computations 
disks are assumed to be slim, which means that the radial advection is taken into account while computing effective temperature of the disk. 
Aims. We attempt to reconstruct continuum and line spectra of X-ray binaries in soft state, i.e. dominated by the disk component of multitem- 
perature shape. We follow how the iron line complex depends on the external irradiation, an accretion rate and a black hole spin. 
Methods. Full radiative transfer is solved including effects of Compton scattering, free-free and all important bound-free transitions of 10 main 
elements. Moreover, we include here the fundamental series of iron lines from helium-like and hydrogen-like ions, and fluorescent K a and 
lines from low ionized iron. 

Results. We consider two cases: non-rotating black hole, and black hole rotating with almost maximum spin a=0.98, and obtain spectra for 
five accretion disks from hard X-rays to the infrared. In non irradiated disks, resonance lines from He-like and H-like iron appear mostly in 
absorption. Such disk spectra exhibit limb-darkening in the whole energy range. External irradiation causes that iron resonance lines appear 
in emission. Furthermore, depending on disk effective temperature, fluorescent iron K B and Kp lines are present in disk emitting spectra. All 
models with irradiation exhibit limb-brightening in their X-ray reflected continua. 

Conclusions. We show, that the disk around stellar black hole itself is hot enought to produce strong absorption resonance lines of iron. 
Emission lines can be observed only if heating by external X-rays dominates over thermal processess in hot disk atmosphere. Irradiated disks 
are usually brighter in X-ray continuum when seen edge on, and fainter when seen face on. 
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1. Introduction 

1.1. General 

The newest X-ray satellite spectra of galactic black hole bina- 
ries (GBHBs) clearly show the existence of UV/ X-ray bump 



peak ing at about 1 keV in their high/soft states (IDavis et al 
2006). It is widely believed, that this emission originates in 
an accretion disk ci rculating around th e central black hol e 
of stellar-like mass dMitsuda et al.l Il984t iDotani et al] [r997). 



Since individual rings of the accretion disk exhibit various 
effective temperatures, the radiation spectrum integrated over 
radii is responsible for multi-blackbody shape of the observed 
spectrum. Similar optical/UV bumps were observed in the 
spectr a of active galactic nuclei (AGN) dShields 1978; Malkan 
19831) . However, accretion disks in AGN exhibit relatively low 



effective temperatures, since the mass of the central black hole 
is of the order of 10 8 M o , and the disk temperature is inversely 



proportional to the mass with the power of 0.25. Multicolor 
disk model (MCD) still does not expla in some spectral details 
jMerloni et alfeOOfi Ikubota et al1l201fjh but up to now there is 
no better explanation for the observed continuum. 

In addition, in both cases of GBHBs and AGN, substan- 
tial fraction of the bolometric luminosity is emitted in hard 
X-ray tail extanding up to several tens or even hundreds keV 
dGierlinski et al.lll997t iNandra & Poundslll994l) . The shape of 
many X-ray tails shows the evidence of interaction between ra- 
diati on emitted near the central object and the nearby accreting 
gas dPounds et al . 1990). The most prominent evidence of such 
an interaction is the detection of fluorescent iron K a line at 6.4 



keV in different tvpes of accretine sources (Tanaka et al. 


1995 


Iwasawa et al. 


1999: Miniutti et al. 2004: Miller et alJ 


2004 


Di Salvo et al. 


2005 


). 
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Shape of an iron line profile carries important information 
about the physics of the material around the compact object. 
Close to the black hole, the emission line profile is relativisti- 
cally broadened and skewed, which together with a circular gas 



2 



Rozanska et al.: Iron lines in model disk spectra of GBHBs 



movement in the di sk makes the final profile asymmetric with 
a strong red wing dFabian et alJll989t iRevnolds & Begelmanl 



1997; Fabian et al 



2002). But, there is a one example of 
broad line profile det ected in cataclysmic variable, GK Per 
dTitarchuk et al.l 120091) . where relativistic effects cannot work 
at all, suggesting that different mechanizm can be responsi- 
ble for the line broadening. Also in so me cases, the detected 



iron line profile has narrow componen t (Watanab e et al. 2003 



Reeves et al 2004; lYaqoob et al. 120071) . indicating that interac 
tion between hard X-rays and accreting gas occurs farther away 
from the black hole. 

In all models used to fit X-ray data, iron line profile is as- 
sumed to be gaussian with eventual kinematic broadening. In 
this paper we present detailed computations of iron line pro- 
files in the case of GBHBs, taking into account thermal, natu- 
ral, pressure and Compton broadening of lines when radiation 
passes through the atmosphere. Suc h calculations were already 
done in case of the disk in AGN dRozanska & Ma dei 2008, 
hereafter RM08), showing profiles of fluorescent K a , Kg lines 
and their Compton shoulders. But accretion disks in GBHBs 
are hot enough for creation of resonance iron lines due to ther- 
mal processes in the atmosphere. The main goal of this paper 
is to show how strong will be full iron line complex (6.3-7.1 
keV) from hot illuminated disk atmospheres in GBHBs. 



1.2. Global disk models 

Numerous theoretical models of the accretion disk spectra were 
recently developed, still a lot of effort reminds to be done to 
achieve undisputable solution. One of the obvious reason for 

such situation is that each mechanism responsible for the ac- 

1 II ^1* 

cretion makes disk matter very turbulent (Ohsuga 2007), and 
no code at present can compute radiative transfer through mov- 
ing gas, no matter if this mechanism is driven by viscosity 
( Shak ura & Sunvaevlll973b . or by magnetorotational instabil- 
itv dHawlev et al.ll200ll) 

Second serious problem is that depending on geometrical 
and physical assumptions we have several global models of ac 



Shakura & Sunvaev 



cretio n disks. The earliest was the standard 

(119731) . SS disk, which is geometrically thin, and opti- 
cally thick, with all dissipated energy converted into ra- 
diati on. Relativistic co r rectio ns for such disk was derived 
by iNovikov & Thornd d 19731) (hereafter NT disk). But in 
other models some thermal energy or mass in the disk 
can be carried into black hole by radial advection. Such 
a process generates two branches of disk solutions: op- 
tically thin ADAFs (advection dominated accretion flows) 
(Nar avan & Yill 1995b. and op tically thick, slim disks developed 
bv lAbramo wicz et al. dl988h. Models of disk spectra presented 
recently by Davis et al.1 d2006h and storred in xspec library as 
fits templates, were computed for NT disk. 

We present, for the first time, sample spectra o f optically 



thick, slim accretion disks ( Abramowi cz et al.lll9 88). We trace 



the behaviour of hot and cold (fluorescent) iron line profiles as 
the function of distance from the black hole, accretion rate, and 
aspect angle. 



1.3. Modeling theoretical disk spectra 

Theoretical disk spectrum can be obtained on various levels 
of sofistication, from the simple black body emission to the 
detailed numerical simulation of the radiation field. In both 
cases locally emitted spectrum has to be integrated over the 
disk surface. The first approach is quite simple as soon as we 
kno w the radial distribution o f effective temperatures in the 
disk (Sha kura & Sunvaevll973l) . The second approach depends 
on the assumptions made in calculations of vertical structure 
and radiative transfer in the disk atmosphere. The most impor- 
tant aspects which have to be taken into account are: treatment 
of energy-dependent opacities of elements, free-free processes 
together with Compton scattering, vertical stratification of tem- 
perature, density and degrees of gas ionization, illumination 
by an external X-ray source and eventual fluorescent emission, 
vertical energy dissipation by accretion, eventual influence of 
vertical height dependent gravity, which operates in the disk. 
Close to the black hole relativistic corrections on the vertical 
disk structure and emitted spectrum should be taken into ac- 
count. 

In the pioneering work done by lLaor & Netzei ( 1989b 
radiative transfer was solved in the Eddington approx- 
imation, neglecting Com p ton s cattering. On the other 



hand, Shimura & Takaharal (1993) considered fully ionized 
hydrogen-helium di sks assuming only bremstrahlung and 
Compton scattering. iDoerrer et al.l (11996b has taken into ac- 
count some absorption but only that of pure hydrogen. None 
of those papers have treated external illumination. 

On the other hand, there were studies of reprocessing of the 
external hard X-rays by a cold matter, neglecting or simplyfy- 
ing the treatment of hydrostatic equillibrium or detailed source 
function calculations. Reprocessed radiation, created mostly 
by Compton scattering process, coul d be succesfully recon 



structed by Monte Carlo simu lations dGeorge & Fabian 1991 



Magdziarz & Zdzi arskJl 1995b . with true absorption taken into 
account in a very schematic way. 

Reprocessing of radiation by Co mpton scattering coul d 
be described by Kompaneets equation IRoss & Fabianl dl993f) ; 

For many years it was the most advanced 



Merloni et 



;d by_l 
aDdlQ 



approach to compute the iron line models in the reflected X- 
ray spectrum. But those authors assume that reflection occurs 
on constant density slab with optical depth of the order of 3, 
while thermalization of hard X-rays may occur deeper inside an 
atmosphere. Hydrostatic equillibrium was a serious comp lica- 
tion of the problem (as discussed by Rozanska et al . 2002, and 
references therein). Taking into account the constraint of hy- 
drostatic equillibrium we can precisely calc ulate strong density 
and i onization gradient in the atmosphere ( Mad ei & Rozanska 
2004b . Even if resulting X-ray spectra from constant density 



slab do not differ significantly in the spectral domain where 
iron line appears ( Ross & Fabianl l2007). we point out here that 
they do not reproduce properly thermal optical UV bump, and 
the strenght of soft X-ray lines, the latter depending on ioniza- 
tion vertical profiles. 

Very sofisticated calculations of non-illuminated accretion 



disk sp ectra were done with tlusty code developed bv lHubenv 
(119900 . They assume disk atmosphere to be in hydrostatic and 
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radiative equillibrium, and solve t he radiative transfer equa - 
tion with non-LTE equation of state ( Hubenv & Hubenvll997 ), 
heavy element abundances an d Compton scatter ing treated 
with the Kompaneets equation (IHubenv et al ] l200ll) . However, 
those authors do not take into account effects of the external 
X-ray illumination, so they cannot precisely compute iron line 
profiles, Nevertheless, the level of sofistication of our radia- 
tive transfer calcul ations is comparable to that presented by 



Hubeny et al. (2001) 



1.4. The aim of our paper 

In this paper, we present precise calculations of the radiative 
transfer through the atmosphere of a slim disk around a stellar 
mass black hole. The global slim disk model ((Sadowski et al 



2010) was computed for several values of the accretion rate and 
for two values of the black hole spin: and 0.98. 



We use here the code developed by Madei (1991) for 
hot stellar atmospheres and then adopted to atmospheres il 



luminated by external hard X-ray source ( Madei & Rozanska 



2000) 



Rozanska & Madei 



and to accretion disk geometry 
200 11) . Adopting atmospheric computations for an accretion 
disk, we solve the vertical structure and outgoing spectra for 
several neighbouring rings being in hydrostatic and radiative 
equilibrium. The final spectrum presented in specific intensity 
scale for different aspect angles is integrated over radii and pre- 
sented at the source frame. No kinematic special relativity ef- 
fects are included in our model. 

Our radiative transfer equation includes effects of multiple 
Compton scattering of radiation on free electrons in relativistic 
thermal motion, and rich set of bound-free and free-free opaci- 
ties. Ionization populations are computed assuming LTE equa- 
tion of state for ideal gas. The method allows for a large rela- 
tive photon-electron energy exchange at the time of Compton 
scattering and, therefore, are able to reconstruct Compton scat- 
tering of photons with energy approaching or e ven exceeding 
the electron rest mass dMadei & Rozans ka 2004). 

We compute profiles of thermal and fluorescent iron lines 
originating from the accretion disk atmosphere irradiated by 
hard X-rays with power-law spectral distribution. Taking into 
account emissivity of fluorescent K ff and lines of low ion- 
ized iron in full NLTE, we are able to estimate the relative im- 
portance of thermal and nonthermal (fluorescent) iron lines. 

In the next section we present the global structure of a slim 
accretion disk model model and principal equations used in our 
radiative transfer computations. Results are presented in Sec. [3] 
and main conclusions are given in Summary section. 

2. Method and model assumptions 

All our models are calculated assuming black hole of the mass 
10M o . Throught all paper we express radial distances, R, in 
units of the Schwarzsild radius defined as rs c hw = 2GM/c 2 . 

In the first stage of calculations we solve global two- 
dimensio nal, hydrodynamical, re lativistic slim disk model pre- 
sented in Sadowski et al. I d2010h . Contrary to the classical SS 
and NT disk models it does not assume radiative equilibrium, 
i.e. exact balance between the dissipated energy and the cooling 
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Fig. 1. Radial dependence of vertical gravity (upper panel) and 
effective temperature (lower panel) in our slim disk models for 
three accretion rates: m =0.01, 0.1 and 1, denoted as black, 
blue and red respectiverly. Two dashed lines (blue and black) 
represent models for wich we assume spin of a black hole 
equaled 0.98. Thin horizontal lines mark the position of hori- 
zont and ISCO for rotating (dashed line) and non-rotating black 
hole (solid line). Points mark annuli for which detiled radiative 
transfer computations were performed. 



rate at a given radius. The slim disks allow for advection of heat 
with inflowing matter which becomes important for moderate 
and high accretion rates (th > 0.3M£^Q) and which affects the 
disk appearence e.g., its effective temperature and pho tosphere 
profiles. The model adopted here (Sad owski et al.ll2010l) solves 
the slim disk equations in two dimensions (i.e. assuming only 
the axial symmetry). The vertical structure is not assumed ad 
hoc but is solved consistently. We assume that the dissipation 
rate is proportional to the local pressure (i.e. we adopt the a 
prescription with a = 0.01). The energy transport in the verti- 
cal direction is treated under the diffusive approximation with 
convection described by the mixing-length theory. The fraction 
of the dissipated heat which is advected, is independent on the 
vertical coordinate by assumption. 

We computed 3+2 models of accretion disks for the accre- 
tion rates m =0.01, 0.1 and 1 in units of Eddington accretion 
rate. For all values of the acretion rate we consider nonrotating 
black hole a = 0, while for two lowest values of m models with 
a = 0.98 were also computed. 

Global models considered in this paper are presented in 
Fig Q] Upper panel represents radial distribution of vertical 
gravity on the disk surface, and lower panel - radial distribu- 



1 Afeid = 16LEdd/c 2 is the critical accretion rate that for a 
disk around a non-rotating BH approximately corresponds to the 
Eddington luminosity, L EM . 
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tion of disk effective temperature. Both quantities are input pa- 
rameters in our radiative transfer calculations. For each global 
model we divide disk into rings centered at the distances shown 
in FigQ~|by rectangle points. For those distances detailed radia- 
tive calculations were performed. 

As the black hole rotates the innermost stable circular or- 
bit (ISCO) and event horizont moves closer to the black hole. 
Therefore, in case of a rotating black hole we can receive also 
radiation from below R=3 (long dashed lines in Fig [TJ. The 
latter value defines the location of ISCO in nonrotating case. 

In all cases, the disk matter is assumed to have solar-like 
chemical composition. Number abundances of elements rela- 
tive to hydrogen, (N e i em /N H ) equal to: 1.00 (for H), 9.54 x 10~ 2 
(He), 4.72 x 10~ 4 (C), 9.65 x 10~ 5 (N), 8.55 x 10~ 4 (O), 
3.84x 10~ 5 (Ne), 4.17 x 10~ 5 (Mg), 4.94x 10" 5 (Si), 1 .64x 10~ 5 
(S), 6.58 x 10~ 5 (Fe). 

2.1. Radiative transfer calculations 

Detailed description of the present code and the corresponding 
equations were given in paper RM08, and here we remind only 
few of them. 

The equation of transfer for the specific intensity I v at fre- 
quency v is solved in plane-parallel geometry on the monochro- 
matic optical depth t v : 

T j y T C / 1 \ 

fl—=I V =Iy-Sy, (1) 

ClTy Ky + CTy 

where S v is the frequency dependent source function. In this 
paper we use the LTE (local thermodynamic equilibrium) ab- 
sorption K v , whereas coefficients of emission j v and scattering 
(T v include non-LTE terms. Emission coefficient j Y is the sum 
of three terms, j v = jf + j v c + j{! , which represent thermal emis- 
sion, Compton scattering emission and the emission in iron flu- 
orescent lines, respectively. 

Coefficient of thermal emission j v in LTE is proportional to 
the Planck function jf = k v B v . The coefficient of true absorp- 
tion K Y , is the sum of bound-free absorption from numerous lev- 
els of atoms and ions for all elements, plus free-free absorption 
from all ions. We also included absorption of 4 lowest lines of 
fundamental series of helium-like iron and of similar 4 lowest 
lines of hydrogen like iron, all formed in LTE by the assump- 
tion. 

The external intensity from the point-like lamp is emitted in 
the form of power-law with spectral index ax, and exponential 
cut-off limits v,„,„ and v max : 



(2) 



(3) 



I ^ = Av -ax expf-— ) expLM], 

\ Vmax) V V I 

and normalized to the luminosity of the source, Lx'. 

A - hL 

4*7 If ex P(-;±) «p(-&*)dv] 

The distance from an X-ray source depends on the ring radius 
r, and the lamp height hi, in the casual relation r 2 = h 2 + r 2 . 
The luminosity, spectral index, and cut-off limits of irradiat- 
ing spectrum are free parameters of our model, all of them de- 
scribed Sec. [3] 



Expressions for the Compton scattering terms and emis- 
sion coefficient in an irra diated atmosphere were derived in 
Madei & Rozanska (I2004I) and are valid in a disk atmosphere 
without any changes 



Jv 



CTy 



+ (Ty 



<S>x{v,V)dV 
o 

(Jy + Uy) <D 2 (V, V) dV . 



(4) 



In the above equation variable U Y denotes the angle- 
averaged intensity of the external irradiation (RM08). 

Compton scatterin g cross sectio ns were computed fol- 
lowing the paper by iGuilberj (11981I) . Functions <l>i and <l>2 
are properly weighted angle-averaged Compton redistribution 
functions for photons both incoming or ou tgoing of frequency 
v afte r scattering in thermal electron gas (Madei & Rozanska 
2004h . 

Fluorescence of low-ionized iron gas was approximated by 
two emission lines, K a and Kg. Therefore, 



f v l = E{ l <p a v +E f B \ 



(5) 



fi ft 

where E a and Ei denote the integrated intensity of K a and Kp 
emission lines, respectively. Frequency dependent variables ip" 
and ijl define profiles of fluorescent lines, both normalized to 
unity. 

Energy emitted in K a and Kp lines was derived by absorp- 
tion of hard continuum X-rays from the radiation field penetrat- 
ing the disk atmosphere. Hard X-ray photon interacting with 
neutral or low-ionized iron most probably ionize and remove 
electron from the innermost K shell, and then remaining hole 
is filled by another electron falling from L shell (K a ) or M shell 
(Kp transition). The integrated emissivity for Fe K a fluorescent 
line is given by: 



Y x hvi 



f 



hv 



■(Jy + Uy)dV, 



(6) 



where we set Y = 0.34 (IKrause & Oliverlll979l) . Variable < on 
is the bound-free absorption coefficient for ionization from K 

shell of iron counted for 1 atom. Similar expression holds for 

fi 

the integrated intensity El . 

The fluorescence yield, Y, defines the fraction of the energy 
of hard X-rays absorbed by iron atoms which was reemitted as 
photons in K a line. The remaining energy of absorbed X-rays, 
1 — Y, was spent for the release of Auger electrons. 

Iron K a fluorescent line is a doublet line, and such a struc- 
ture was reproduced by our code. The doublet structure is 
caused by the fact that L shell has three sub-shells depending 
on the value of the spin and orbital angular momentum, and 
only two transitions are allowed by selection rules. We set the 
central energies for K ai and K a , lines to 6.404 keV and 6.391 
keV, respectively. Natural widths (FWHM) of both lines, 2.7 
eV (K a ,) and 3.3 eV (K an ) were taken from IKrause & Oliver 
(11979b . The Kp line was approximated by a singlet line at the 
energy 7.057 keV. We set the natural width of the line to 2.5 eV 
arbitrarily. 

Integrated emissivity E{, was divided in the proportion 2: 1 
between K ai and K a , components of the doublet, since the 
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lo g( T «) lo g( T eJ l0 g( T es) 

Fig. 2. Vertical temperature structure for various rings in three non-illuminated disk models: Models differ by accretion rate m 
and dimentionless spin of the black hole a. Curves are labeled by the distance from the black hole, R. In all cases we observe 
temperature rise in outermost layers of atmosphere due to Compton scattering (see text for details). 



number of electrons on L shell with total spin equaled 1 is 
twice higher than the numb er of electrons with total spin equal 
0. Following iBaskol (11978b we quite arbitrarily assumed that 
the integrated intensity E^ = 0.1 E^. Opacity profiles of all 
three lines were set to Voigt functions with depth-dependent 
parameters describing natural and thermal broadening. 

On each discrete ring we performed radiative transfer cal- 
culations assuming that the disk atmosphere remains in hydro- 
static and radiative equillibrium (Eq. 11 and 12 in RM08). No 
turbulences and convection were taken into account. 

Full set of model atmospher e equations is solved b y 
the method of partial linearization (IMadei & Rozans ka 2004). 
Equation of hydrostatic equillibrium was excluded from the 
linearization scheme. The equation of transfer, Eq. [T]is solved 
by the Feautrier meth od and the variable Eddington factors 
method dMihalasll 1978b . 



Table 1. Description of lines in Figs. |4]-[9] 



3. Results 

Our code allows us to compute the structure of disk atmo- 
spheres over very large range of electron scattering optical 
depth starting from r es - 10~ 8 up to r es = 10 5 . Furthermore, 
we were able to reproduce the overall continuum spectrum 
from deep infrared of 0.4 eV up to hard X-rays of 400 keV. We 
solve the radiative transfer problem on the grid of 175 optical 
depth and 2200 photon energies simultaniously. 

All spectra are presented as energy dependent outgoing 
specific intensities, Ie, which are suitable for disk geometry. 
We reject presentation of monochromatic fluxes since they are 
relevant only to geometry of a spherical star. Our intensities 
are shown for 8 cosines of viewing angles, which represent an- 
gles of the Legendre quadrature. Exact values of those angles, 
their cosines and the type of lines used in all figures presenting 
disk spectra are given in TableQ] In further discussion we draw 



type of a line 


cos(/) 


i 


solid black 


0.9801 


11.4° 


short-long dashed 


0.8983 


26.1° 


long dashed dotted 


0.7628 


40.3° 


short dashed dotted 


0.5917 


53.7° 


dotted 


0.4083 


65.9° 


short dashed 


0.2372 


76.3° 


long dashed 


0.1017 


84.2° 


solid red 


0.0199 


88.9° 



attention of the reader to the extreme angles: solid black line 
represents almost vertical direction (face-on aspect), whereas 
the solid red line represents almost horizontal direction (edge- 
on aspect). 

We consider only one choice of irradiating power-law with 
ax = 0.7 and luminosity Lx = 10 35 erg s , ranging from 
hv m i„ — 0.1 up to hvmax = 100 keV. Position of the X-ray lamp 
is fixed arbitrarily on the height hi = 5rs c hw above the disk at 
the radius R — 2 from the black hole. 



3.1. Temperature structure 

Radial distribution of effective temperatures in all global mod- 
els has characteristic shape (see Fig. Q] For non rotating black 
hole the maximum value of the effective temperature occurs 
on R « 5. Closer to the black hole, there is a slight temper- 
ature drop, which is stronger for the lower accretion rate. For 
m = 0.01, the effective temperature at R=3 has the same value 
as on R about 60, while the surface gravity is sligtly higher. 
This matters in case of irradiated atmospheres, where the flu- 
orescent line can be created due to illumination of innermost 
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Fig. 3. Comparison of the temperature structure in three disk models before (black line) and after (red line) illumination by 
external radiation in the form of power-law. Results are shown for various distances from the black hole. It is clearly seen that 
the external illumination affects only outermost layers of the disk atmosphere. The most prominent effects of irradiation occur 
at distances where the effective temperature of an atmosphere itself is relatively low. It happens always when the accretion rate 
is low, see left panel. For high accretion rate, high spin and close to the black hole, the atmosphere itself is hot and the external 
illumination can slighly cool the outermost layers due to Compton cooling. 



ring with the effective temperature of the order of 7 x 10 5 K 
(see Sec. 13.21 1. Moving farther from the black hole temperature 
decreases, and at R=100 is half order of magnitude lower than 
its maximum value. When the black hole rotates the tempera- 
ture of innermost rings increases reaching 10 7 K, and we can 
expect high ionization degree of iron in such conditions. 

Vertical temperature profiles computed for individual rings 
are presented for various global models and distances at Fig. [2] 
and[3] All profiles of the disk atmospheres exhibit inverse tem- 
perature rise in the uppermost atmospheric layers. In the case of 
non-illuminated disks, presented in Fig. [2] the inverse tempera- 
ture rise is due to Compton heating by hot X-ray photons from 
below the photosphere. This result is in agreeme nt with models 
of hot neutron star atmospheres com puted by iMaiczvna et al. 

Sulei manov & Wernerl (I2007|). and with atmosphere of 
accretion disks presented bv lHubenv et alJ d200ll) . Note that in 
Fig. |2] we present three hottest models where the surface tem- 
perature reaches 2.2 x 10 7 K. In such conditions, iron is almost 
fully ionized due to heating by radiation generated already via 
accretion process. 

In disks with external illumination the heating effect is en- 
hanced (see red lines in Fig. [3]) due to Compton scattering of 
hot external X-ray photons. The amount of this heating depends 
on the disk effective temperature determined by the accretion 
rate and the black hole spin. For the lowest accretion rate and 
non rotating black hole (left panel) the effect of irradiation is 
seen in all rings, since effective temperatures for this model are 
the lowest. When the black hole rotates (middle panel) X-rays 
heat up mostly outer rings. The same is observed for higher ac- 
cretion rate and non rotating black hole (Fig.[3]right panel). We 
do not show temperature run for two hottest disks, since then 



the external irradiation does not cause any significant heating 
effect. 

Our results include the influence of both Compton heating 
sources: by hot photons from the photosphere, and from the ex- 
ternal X-ray source. Because of that, surface temperature struc- 
ture depends on the shape and the normalization of the external 
irradiation. In case when the atmosphere is hot by itself, addi- 
tional soft X-ray photons can cool down the hottest layers due 
to Compton down-scattering. This happens for rings located 
close to the black hole, for the model with m = 1 and a = 0, 
where the temperature of the uppermost irradiated surface lay- 
ers is lower than for non-irradiated atmosphere (Fig. [3] right 
panel). Compton down-cattering is mostly caused by photons 
with energies lower than gas temperature, so we can expect that 
hot disk atmospheres in galactic black hole binaries can be cool 
down due to irradiation by soft X-rays. We plan to study this 
effect in a forthcoming paper. 

3.2. Spectra from individual annuli 

Intensity spectra emitted by individual rings of non irradiated 
accretion disks are presented in Fig. [4] Each panel shows evo- 
lution of resonance lines from He-like and H-like iron for four 
models with the lowest effective temperatures. Line profiles 
change with distance from the black hole, which is marked on 
the right vertical axis. In some rings the line at 6.7 keV ex- 
hibits reverse emission in the line core which (for LTE) is the 
signature of the temperature inversion. 

For the disk with the lowest accretion rate and zero spin, 
only absorption line from He-like iron is present originating 
in rings of the maximum effective temperature. Spectra from 
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m=0.01; a=0 fn=0.01; a = 0.98 m=0.1; a = fn=l; a = 




6.7 6.8 6.9 6.7 6.6 6.9 6.7 6.8 6.9 6.7 6.8 6.9 

Energy (keV) Energy (keV) Energy (keV) Energy (keV) 



Fig. 4. Outgoing intensity spectra from individual rings in four non illuminated disk models at different accretion rates and black 
hole spins. Spectra are shown in the energy range corresponding to the resonanse lines from hydrogen and helium like iron. 
Spectra are shifted along vertical axis in order to demonstrate evolution of lines with distance from the black hole (marked on 
the right vertical axis). The limb-darkening effect, where intensity from "face-on" disk (black solid line) is higher than intensity 
from "edge-on" disk (red solid line), is observed in all cases, except in some line cores. 



m = 0.01; a = 0; a = 0.7 fn = 0.01; a=0.98; a = 0.7 fn = 0.1; a = 0; a = 0.7 




6.4 6.6 6.8 7 6.4 6.6 6.8 7 6.4 6.6 6.8 7 

Energy (keV) Energy (keV) Energy (keV) 

Fig. 5. Outgoing intensity X-ray spectra for sample of irradiated accretion disks. Spectra are shifted along vertical axis in order 
to demonstrate the evolution of iron line complex with distance from the black hole (marked on the right vertical axis). 



two rings located near the black hole, and from all rings above 
R = 20 are featureless, since temperature is low and there is not 
enought He-like and iron ions. Situation changes for the same 
accretion rate when the black hole rotates. Iron line at 6.7 keV 
is stronger and the line at 6.9 keV becomes visible. For m — 1 
and a — both lines are clearly seen for most of radii R. There 
exists energy dependent limb-darkening in all non irradiated 
models. At some rings, differences between intensities emitted 



at various angles are strong, while sometimes there are not seen 
at all. 

Spectra of irradiated disks are presented in Fig.|5]for the en- 
ergy range 6.3-7. 1 keV. For most rings both resonance lines of 
iron appear in emission. This effect, caused by external illumi- 
nation is visible only for those rings, when there is temperature 
change due to irradiation (see also Fig. [3]). For the model pre- 
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Fig. 6. Broad-band intensity spectra integrated over disk sur- 
face for all non-illuminated disk models for eight viewing an- 
gles. The shape of continuum is affected by several spectral 
edges and iron line complex around 6.9 keV. Depending on the 
accretion rate and spin of the black hole we can follow how 
limb darkening changes with energy. 



sented in the middle panel for R < 6, both resonant lines appear 
in absorption. 

Rings with relatively lower effective temperature, usually 
located farther from the black hole show cold fluorescent iron 
lines, K„ at 6.4 keV and Kg at 7.05 keV. We point out here that 
those lines can be visible in the whole iron line complex, when 
integrated over disk surface (see the following subsection). 

External irradiation causes limb-brightening in intensity 
spectra emitted from numerous rings. In the energy domain 
of Fig. [5] continuum limb-brightening apparently is correlated 
with the reverse emission in both resonance lines of iron. For 
the model presented at the middle panel, lines emitted from in- 
nermost rings are in absorption and simultaneously we observe 
usual limb-darkening. 



3.3. Disk-integrated spectra 

At the last step of computations we integrate spectra emitted 
from individual rings over disk surface according to the usual 
formula: 



l'°'{ f i) = 27r)\l ( "\ij)R„AR„ 



/ j 1 E 
n=l 



(7) 



where n is a ring number, // is the cosine of viewing angle, and 
AR n is the width of the n-th ring. 



Fig. 7. Profiles of resonance lines from helium-like and hy- 
drogen like-iron, integrated over disk surface for all non- 
illuminated disk models. No kinematic either general relativity 
effects were taken into account. 



Fig. |6]presents continuum spectra integrated over disk sur- 
face for all non irradiated models. Left hand part of each panel 
displays featureless UV, optical and near infrared spectrum of 
disk whereas right hand part reaches hard X-rays and iron line 
region. All spectra show limb-darkening over the whole range 
of energies. This is clearly seen in Fig. |7j intensity edge on 
(solid red line) is much lower than face on intensity (solid black 
line). 

Resonance iron lines appear always in absorption for non 
irradiated models (Fig. |7J. They are very deep, i.e central in- 
tensity in the line can drop almost three orders of magnitude 
(panel a) as compared to the continuum level. Our code takes 
into account pressure broadening of all resonance lines, and 
particulary Lyo- line of He-like iron can be very wide, as seen in 
panels c and e. Such prominen t lines certainl y can be detected 
in present X-ray observations (IKubota et al]l2007l) . Accretion 
disks with highest effective temperatrue do not show strong 
iron lines in absorption, and for the hotest model iron is fully 
ionized, as seen in panel d. 

Continuum spectra of the same disk models with exter- 
nal irradiation are displayed in Fig. [8] All spectra exhibit 
Compton reflection bump in hard X-rays. The relative posi- 
tion of Compton reflection in comparison with UV/X-ray disk 
bump depends on accretion rate and black hole spin. External 
irradiation causes shift to limb-brightening in the whole region 
of Compton reflection bump. For low accretion rate and non ro- 
tating black hole limb-brightening appears also in visible and 
near infrared spectral region (panels a and c in Fig . [8} . 
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Fig. 8. Broad-band intensity spectra integrated over disk sur- 
face for all irradiated disk models for eight viewing angles. 
The shape of continuum is extended toward hard X-rays due 
to both: thermal emission of hot atmospheric zone and the ra- 
diation from the external source reflected by Compton scatter- 
ing. It is evident that continuum intensity spectra exhibit limb- 
brightening in X-ray domain. 



Iron line complex for illuminated accretion disk models is 
presented in Fig. [9] Structure of the iron line complex is rather 
complicated and its properties, i.e. line absorption/emission 
and limb darkening/brightening, do not depend in the straight- 
forward way on the accretion rate and the black hole spin. 

Lines from He-like and H-like iron appear in emission only 
for low accretion rates. For the accretion rate m = 1 those lines 
are in absorption. Higher level of irradiation may change this 
into emission, but this depends on the shape of the illuminat- 
ing radiation. Increasing of soft X-ray photon number leeds to 
the efficient Compton down-scattering, which effectively cools 
down already hot atmosphere in galactic black hole binaries 
(see Sec. 13. 11 1. 

Resonance lines in emission presented in panel b are very 
wide and do not show any limb brightening. We suppose that 
this wide profile is caused by Compton scattering in hot atmo- 
sphere. As expected, model presented in panel d is featureless. 

The important result of our paper is that in models with 
lower effective temperature, cold fluorescent K a and lines 
are emitted. The strenght of these lines is large enough to con- 
tribute to the final iron line complex. In our code we assume 
that the cold fluorescent K a line has fixed energy 6.4 keV. But 
this energy depends on detailed ionization structure of iron 
atoms, and we shall discuss this problem in a forthcoming pa- 
per. 




6.4 6.6 6.8 

Energy (keV) 



Fig. 9. Iron line complex in the energy range 6.3-7.2 keV, in- 
tegrated over disk surface for all illuminated disk models. In 
three upper panels, Lyman alfa lines of helium-like iron (6.7 
keV) and of hydrogen-like iron (6.9 keV) are in emission, 
which is caused by the layer of temperature invertion in the 
uppermost atmosphere. Only for the hottest disk models (two 
bottom panels) both lines are in absorption, since the layer of 
temperature inversion is thin. Fluorescent emission lines, K„ at 
6.4 keV and Kp at 7.05 keV, can be seen only in colder disk 
models. 

4. Summary 

In this paper we present model atmospheres and intensity spec- 
tra for a set of accretion slim disks around black holes in galac- 
tic black hole binaries. Model atmospheres are computed for 
various accretion rates and spin of the black hole. Our com- 
putations include effects of Compton scattering and external 
irradiation by an X-ray lamp. Spectra are obtained for a large 
number of individual rings, depending on the particular model, 
the average number is 14. Outgoing intensity spectra were com- 
puted for wide energy range from infrared to hard X-rays, up 
to 400 keV. Note, that results presented in our previous paper 
(RM08) were calculated for SS disk global models, whereas 
the present paper deals with slim disks with advection. 

We demonstrate that the external irradiation developes huge 
temperature rise in the outermost layers of disk atmospheres. 
The effect is most prominent in cases when the atmosphere 
has lower effective temperature (in the outer and the innermost 
rings of our slim disk models). In those colder rings, strong 
Kq, and Kjg fluorescent lines of iron are emitted. Results of 
the paper demonstrate also that the external irradiation causes 
a change of limb-darkening of non irradiated disk to limb- 
brightening particularly in hard X-ray region. This implies that 
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accretion disks seen edge on are brighter in X-rays than disks 
seen face on. 

Our models clearly show that for high accretion rates it is 
rather difficult to obtain strong hot iron resonance lines in emis- 
sion. This is due to the fact, that the disk atmospheres in galac- 
tic black hole binaries are generaly very hot. In such a situation 
iron can be fully ionized and produces no line features. For low 
accretion rates, resonance iron lines appear in emission, when 
the heating by comptonization of external photons causes huge 
temperature rise. 

Our results differ from previously computed spectra 
from constant density slabs of 10 Thomson optical depth 
(Ros s & Fabianl l2007t iGarcfa & Kallmanl 120101) . This papers 
did not take into account the possibility of the resonance iron 
lines formation in hot atmospheres of the accretion disk around 
a stellar black hole. In both cases iron lines appear only in emis- 
sion. 

Results of this paper show that both resonance iron X-ray 
lines can be form in absorption for hot accretion disk atmop- 
sheres in GBHc. This is another possibility to explain obser- 
vations of those ab sorption lines in some blak hole transients 
dKubota et al.l l2007, and references therein). 

Moreover, iron line profiles computed at source must be 
always pressure broadened, and their profiles deviate from the 
gaussian shape. Most of emission and absorption line profiles 
are not gaussian (Fig. [7] and |9). For the same accretion rate, 
iron line profiles differ for two values of spin considered in our 
paper (panel a, b, and c, c in Fig. [9] This should be taken into 
account when fitting of black hole spin to observed the iron 
lines. 

Acknowledgements. This research was supported by the Polish 
Ministry of Science and Higher Education grant No. N N203 5 1 1638. 

References 

Abramowicz, M. A., Czerny, B., Lasota, J. R, & Szuszkiewicz, 

E. 1988, ApJ, 332, 646 
Basko, M. M. 1978, ApJ, 223, 268 

Davis, S. W., Done, C, & Blaes, O. M. 2006, ApJ, 647, 525 
Di Salvo, T, Iaria, R., Mendez, M., et al. 2005, ApJ, 623, L121 
Doen-er, T, Riffert, H., Staubert, R., & Ruder, H. 1996, A&A, 
311,69 

Dotani, T, Inoue, H., Mitsuda, K., et al. 1997, ApJ, 485, L87+ 
Fabian, A. C, Rees, M. J., Stella, L., & White, N. E. 1989, 

MNRAS, 238, 729 
Fabian, A. C, Vaughan, S., Nandra, K., et al. 2002, MNRAS, 

335, LI 

Garcia, J. & Kallman, T. R. 2010, ApJ, 718, 695 
George, I. M. & Fabian, A. C. 1991, MNRAS, 249, 352 
Gierlinski, M., Zdziarski, A. A., Done, C., et al. 1997, 

MNRAS, 288, 958 
Guilbert, P. W. 1981, MNRAS, 197, 451 
Hawley, J. F, Balbus, S. A., & Stone, J. M. 2001, ApJ, 554, 

L49 

Hubeny, I. 1990, ApJ, 351, 632 

Hubeny, I., Blaes, O, Krolik, J. H., & Agol, E. 2001, ApJ, 559, 
680 



Hubeny, I. & Hubeny, V. 1997, ApJ, 484, L37+ 

Iwasawa, K., Fabian, A. C., Young, A. J., Inoue, H, & 

Matsumoto, C. 1999, MNRAS, 306, L19 
Krause, M. O. & Oliver, J. H. 1979, Journal of Physical and 

Chemical Reference Data, 8, 329 
Kubota, A., Done, C., Davis, S. W., et al. 2010, ApJ, 714, 860 
Kubota, A., Dotani, T, Cottam, J., et al. 2007, PASJ, 59, 185 
Laor, A. & Netzer, H. 1989, MNRAS, 238, 897 
Madej, J. 1991, ApJ, 376, 161 
Madej, J. & Rozanska, A. 2000, A&A, 363, 1055 
Madej, J. & Rozanska, A. 2004, MNRAS, 347, 1266 
Magdziarz, P. & Zdziarski, A. A. 1995, MNRAS, 273, 837 
Majczyna, A., Madej, J., Joss, P. C., & Rozanska, A. 2005, 

A&A, 430, 643 
Malkan, M. A. 1983, ApJ, 268, 582 

Merloni, A., Fabian, A. C., & Ross, R. R. 2000, MNRAS, 313, 
193 

Mihalas, D. 1978, Stellar atmospheres /2nd edition/ (San 
Francisco, W. H. Freeman and Co., 1978. 650 p.) 

Miller, J. M., Fabian, A. C, Reynolds, C. S., et al. 2004, ApJ, 
606, LI 31 

Miniutti, G., Fabian, A. C, & Miller, J. M. 2004, MNRAS, 
351,466 

Mitsuda, K., Inoue, H, Koyama, K., et al. 1984, PASJ, 36, 741 
Nandra, K. & Pounds, K. A. 1994, MNRAS, 268, 405 
Narayan, R. & Yi, I. 1995, ApJ, 444, 231 
Novikov, I. D. & Thorne, K. S. 1973, in Black Holes (Les 

Astres Occlus), ed. C. Dewitt & B. S. Dewitt, 343-450 
Ohsuga, K. 2007, ApJ, 659, 205 

Pounds, K. A., Nandra, K., Stewart, G. C, George, I. M., & 

Fabian, A. C. 1990, Nature, 344, 132 
Reeves, J. N., Nandra, K., George, I. M., et al. 2004, ApJ, 602, 

648 

Reynolds, C. S. & Begelman, M. C. 1997, ApJ, 488, 109 
Ross, R. R. & Fabian, A. C. 1993, MNRAS, 261, 74 
Ross, R. R. & Fabian, A. C. 2007, MNRAS, 381, 1697 
Rozanska, A., Dumont, A.-M., Czerny, B., & Collin, S. 2002, 

MNRAS, 332, 799 
Rozanska, A. & Madej, J. 2001, Advances in Space Research, 

28,433 

Rozanska, A. & Madej, J. 2008, MNRAS, 386, 1872 
Sadowski, A., Abramowicz, M., Bursa, M., et al. 2010, ArXiv 
e-prints 

Shakura, N. I. & Sunyaev, R. A. 1973, A&A, 24, 337 
Shields, G. A. 1978, Nature, 272, 706 
Shimura, T. & Takahara, F. 1993, ApJ, 419, 78 
Suleimanov, V. & Werner, K. 2007, A&A, 466, 661 
Tanaka, Y, Nandra, K., Fabian, A. C, et al. 1995, Nature, 375, 
659 

Titarchuk, L., Laurent, P., & Shaposhnikov, N. 2009, ApJ, 700, 
1831 

Watanabe, S., Sako, M., Ishida, M., et al. 2003, ApJ, 597, L37 
Yaqoob, T, Murphy, K. D., Griffiths, R. E., et al. 2007, PASJ, 
59, 283 



